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Qh! Abstract 



We present, in the context of dimensional regularization, a prescription to 
renormalize Feynman diagrams with an arbitrary number of external fermions. This 
prescription, which is based on the original t'Hooft-Veltman proposal to keep exter- 
nal particles in four dimensions, is particularly useful to define the 'renormalization' 
(in the context of effective Lagrangian) of physical four-quark operators without in- 
troducing any evanescent operator. The results obtained for b — > s processes agree 
with those from the so-called naive prescription, but disagree with the ones with 
the introduction of evanescent operators in a renormalization group analysis. We 
also present an explicit two loop calculation of the mixing of the evanescent oper- 
ators with the physical dimension five operators for the same processes. Particular 
attention is paid to the unboundedness nature of such mixing and how a formal 
finite transformation is effected to decouple. The inevitable mass dependence of 
one of these schemes in the literature is pointed out as the cause for the difference 
mentioned. 

PACS numbers(s): ll.10.Jj, 12.38.Bx, 13.40.Hq 
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1. Introduction 

There are at least three independent complete calculations of b — > s processes 
which disagree with each other Refs. [1,2,3]. Although the difference is numerically 
insignificant for b — > S7 (within 2%), it is theoreticaly important to obtain the 
correct Wilson coefficients which contain all the heavy top and W-boson effects. 
The disagreement comes mainly from the introduction of evanescent four quark 
operators in Ref.[2,3]. These were not included in the analysis of Ref.[l], the con- 
tributions of which were subsequently given in Ref.[4], which brought Ref.[l] and 
Ref.[3] into agreement. In this article, we shall attempt to give some partial answer 
to where and when one should introduce these evanescent operators. We shall argue 
that if one is interested in constructing an effective Lagrangian to a particular finite 
order in the strong QCD coupling constant g s , based on direct Feynman diagram 
evaluation of processes, then inclusion or exclusion of these operators is an arbi- 
trary choice. We shall present a scheme for this purpose, in which there is no need 
to include any evanescent operators. It is based on the original t'Hooft-Veltman 
(HV) prescription to keep the external particles in four dimensions. However, if one 
is to perform a renormalization group analysis, then one has to be concerned also 
with the closure of the operator basis under multiplicative renormalization and the 
light mass independence aspect of the Wilson coefficients. Although we have not 
been able to pinpoint, we hope that the following presentation of our concern with 
respect to these last two items may spur further investigation. 

Before we proceed any further into the technical details, we now give a brief 
exposition of the issues, so that an interested reader who may not be a practitioner 
can have an overview. 

A singular development in theoretical physics in the past twenty years or so has 
to do with the successes of formalisms in short distance analysis. In many physical 
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situations, it has been shown that scales far different from the ones under consider- 
ation can be, very loosely speaking, 'integrated out' and their residual effects can be 
parameterized in a series of terms made of effective coefficients and local products 
of low excitations. Particularly in particle physics, let ip be the fields which are 
being prepared and detected experimentally. One is interested here in expanding a 
theory capable of more detailed description at a higher energy scale A into terms 
schematically as 

L full = L low energy + J_ £ Q.Q. + . . ^ 

i 

If L,f ul1 is a renormalizable theory, then L low ener 9y will include all the renormaliz- 
able terms after such an "integrating out' of the heavy particles. Ci are the Wilson 
coefficients, which are functions of log A and perhaps some other variables and Oi 
are made of local functions of ip. We have put a bar over Ci to indicate that their 
factor has been taken out. This is a well-defined construct under the name of 
effective Lagrangian. 

Let us discuss some aspects of the actual construction of an effective La- 
grangian. Specifically, we consider only the case when L^ uU is known and renormal- 
izable. One can perform a loop expansion with it, from which the right hand side of 
Eq.(i) can be extracted. The most systematic way to do this to a particular order is 
by using Zimmermann oversubtraction identities to rearrange the integrand for any 
diagram in which there is at least one heavy internal particle. Terms are organized 
or discarded according to dimensional counting. Since the role of the identities is 
merely to add and to subtract the same terms to a given integrand to allow power 
counting down to parts of it, there are many possibilities. The common tread is an 
expansion in external momenta and internal masses of these subdiagrams relative 
to A. A point which needs to be emphasized, however, is that after the required 
wave function and parameter renormalizations inherent in L^ ul1 , all these parts are 
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finite quantities. If the regulator used for the said renormalizations is dimensional 
continuation, then the limit to go to four dimensions can be taken at this stage 
and each term remains finite. The advantage here is that there are fewer operators 
Oi, which need to be introduced. They are called dimensionally physical. On the 
other hand, one may choose to continue working slightly away from four dimensions 
after each finite order, because the algebraic aspect of the operators may be more 
systematic. This comes about because the mechanics of constructing an effective 
Lagrangian is such that one needs to reinsert the lower order pieces into a diagram 
in order to obtain from it the higher order pieces. In so doing, one may find it 
convenient to introduce extra operators, which are termed evanescent because they 
vanish formally in the four dimension limit. These are some of the common choices 
one is free to make. They all should produce exactly the same S-matrix elements, 
which are dictated by L^ uli , to the accuracy in any inverse power of A demanded 
by us. 

We now turn to the determination of Ci in infinite partial sums. As it is well- 
known, when we include QCD, we should view Ci as functions of g 2 s and g 2 s log A 2 , 
and perhaps of other variables also as we shall see. Because g 2 s log A 2 is not small, 
an improved perturbation expansion is to sum infinite series in leading powers, sub- 
leading powers, etc. of it. The tool for this is the renormalization group analysis. 
Essential to the procedure is the introduction of a dimensionful parameter fj,, which 
in dimensional continuation is associated with the regularization program. As our 
earlier discussion reveals, terms such as J2 i CiOi are finite quantities, upon express- 
ing them in renormalized fields, masses and coupling constants. It should also be 
noted that once a choice is made of the Zimmermann rearrangement, there ensues 
a particular renormalization prescription for O;. Schematically, using the closure 
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property 



% = ZyO] 



bare 



(ii) 



we write 



v-rOi + nfijOj = o, 



(Hi) 



with 



kj > 




and then it follows that 




(v) 



It is obvious that a preferential scheme is one which is mass- independent, such 
that are functions only of # s (/i) and Ci only of g s (n) and log(A///). This allows 
one to perform infinite sums over leading logarithm, subleading logarithm, etc more 
conveniently. 

We shall develop this program in much greater details in the sections to come. 
Thus, the article is organized as follows. In section 2, we review again the general for- 
malism for calculating directly the Wilson coefficients from Feynman diagrams and 
explain how the effective theory is derived from the full theory (Standard Model). 
We also discuss the scheme-dependence of various renormalized quantities in the ef- 
fective theory. In section 3, we present the renormalization group equations (RGE) 
for the Green's functions and the Wilson coefficients and discuss the meaning of the 
leading logarithm approximation (LLA) and the next-to-leading logarithm approx- 
imation (NLLA). In section 4, we present a new prescription to obtain the physical 
four quark operators without introducing any evanescent operators. This scheme 
was implicitly used in Ref.[l] and therefore, if adhered to for (RGE) analysis, all 
the results in the LLA will be the same as those of Ref.[l]. In section 5, we follow 
others in introducing evanescent operators and perform a two loop calculation of 
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the mixing of these operators with the physical dimension five operators. We shall 
discuss the potential difficulty in dealing with unbounded infinite matrices. We 
shall give the similarity transformation which formally decouples these two sets of 
operators at any running scale fx. It is the same transformation which makes the 
matrix elements of these evanescent operators to go to b — > S7, sG vanish at one 
loop. In section 6, we comment on the results and give conclusions. 

2. Effective theory 

Consider the process b — > S7 in the Standard Model (SM). The diagrams 
contributing to this process to lowest order are shown in Fig. 1. Suppose we 
want to calculate the value of each renormalized diagram in the limit m t ~ m w ^> 
p ex t,mb,m s ,m u ,m c to an accuracy C(l/m^), where p ex t represents here the mo- 
memta of the external particles. We use the first diagram, redrawn in Fig. 2, to 
explain how the effective theory is built to this order. We do this by partitioning the 
diagram into several pieces to extract individually the effects, an approach 

which is explained in details in Ref. [5]. Let V denote the unrenormalized diagram 
on the left of Fig. 2 and let us define: 



where p represents the external momenta of the graph (or subgraph) 7, and e = n— 4, 
n being the number of space-time dimensions. The renormalized V is obtained by 
subtracting V at p ext = 0, 



Let 70 be the tree subgraph of V which contains the VF-boson line as illustrated in 



A e) = Pole Part of 7, 



7-W = 7 (P = 0) +Pg-l(P = 0) + • • • + ^P r g-pl(P = 0). 



(1) 




(2) 
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Fig. 3. We can now rewrite Eq. (2) as: 



+ 



+ 



(i-O r<?r 

(r< 2 »-r<°»)r 
(l-r« 2 ')(l-r<»))r. 



(3) 



The first and second terms on the right in Eq. (3) correspond to the first and 
second diagrams on the right in Fig. 2. The last term in Eq. (3) is of order 
l/m^ and is discarded. Eq. (3) is nothing but an algebraic rearrangement, first 
used by Zimmermann [6]. We would like to point out that every term is finite and 
well defined; in particular when computing the quantities in Eq. (3), the Dirac 
algebra can be done with no ambiguity in four dimensions . This rearrangement or 
partitioning of graphs tells us how to compute the heavy effects in the full theory. 
It also generates an effective theory, order by order in perturbation theory. 

In the language of effective theory, the graphs on the right in Fig. 2 are equal 
to the sum of the products 



where the Ci are the Wilson coefficients, and contain all the heavy W-boson and/or 
top quark effects (also called the short distance effects) and the 0\ en are local 
renormalized operators of dimension equal to at most six. The subscript light is 
there to remind us to use light fields (u, d, s, c, 6, G, A) only for the internal as well 
as external lines of the graphs when computing the Green functions with insertions 
of operators. In this case, t$ gives us the tree vertex of Fig. 3, while (1 — t^)t$T 
is the one loop renormalized matrix element of this same vertex, which is the first 
graph on the right hand side of Fig. 2. (t^ — Tp )r denoted by the second graph 
is a part of the effective Lagrangian induced at one loop order, where Ci have been 
evaluated to one loop and r(Oj) have tree values. A general remark is in order. 
Such evaluations of Ci give us the boundary conditions for RGE. The matching 




(4) 
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between the high energy theory and the low energy theory has been automatically 
done. 

We would like to digress for a momemt and discuss the renormalization of the 
Feynman graphs in the full theory and the renormalization of the operators in the 
effective theory. In Eq. (2) we have used momentum subtraction for the graphs in 
the full theory. In Eq. (3) we have also used momentum subtraction to renormalize 
the operators: this is shown in the first term on the right in Eq. (3), where 



prescription of the operators in the effective theory may be taken to be different from 
that of the graphs in the full theory. For example, we may switch to the minimal 
subtraction to define the renormalized operators. This is achieved by rewriting Eq. 
(3) as (neglecting l/m% effects) 



What we learned from this exercise is that one can choose any scheme to renormalize 
the operators. This freedom of choice is also reflected in the dependence of the 
Wilson coefficients: the second term on the right in Eq. (1) and that in Eq. (5) 
depend on the renormalization chosen for the corresponding first term. However, 
the sum CjT\^ ht {p r i en ) is uniquely given once the renormalization of the graphs 
in the full theory is fixed: it equals T ren . One may wonder why we have used 
momentum subtraction for the full theory instead of minimal subtraction. The 
reason is that power counting is easy in momentum subtraction. If we had used 
the MS scheme for the full theory, then some of the C; would have been of order 
m° , and we would not have had explicit decoupling of the W-boson. The effective 
theory obtained would have necessarily required multiple insertions of the associated 




is the counterterm to the bare graph t!$ F. Actually, however, the renormalization 




(5) 
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operators, because they would not have been suppressed by inverse powers of m w . 
In the end, all these multiple insertions could be reabsorbed into the redefinition of 
the fields and parameters in the full theory and one would obtain the same on-shell 
results. On the other hand, using momentum subtraction for the heavy graphs 
insures us that all the Wilson coefficients are of order 1/m^, and therefore we have 
explicit decoupling. The operators are inserted once only. Also the derivation of 
renormalization group equations is simpler when there is explicit decoupling. 

Let us now consider another example which creates " problems" in the effective 
theory. This is illustrated in Fig. 4, where again we want to work to an accuracy 
of I/to^,. Let T be the graph on the left in Fig. 4, and 70 be the tree subgraph 
of T which contains the W-boson line, similar to Fig. 3. V is finite and needs no 
renormalization. The Zimmermann algebraic rearrangement gives 



The first and second terms on the right in Eq. (6) correspond to the first and 
second diagrams on the right in Fig. 4 respectively. The last term in Eq. (6) is of 
order and will be neglected. Here we have again tentatively used momentum 

subtraction to define the renormalized Green functions with 0[ en inserted. In a 
moment we will rewrite Eq. (6) in the MS scheme, but first we would like to make a 
remark. The first and second terms on the right in Eq. (6) are well defined and the 
Dirac algebra can be done in four dimensions safely since each term is finite. This 
means that r ren ((3[ en ) is unambiguously defined, as long as we evaluate the terms 
in the combinations as shown. However, we may want to compute the 'countert- 
erms' ( — V ) separately by themselves, then we encounter in dimensional 




(0) 



(6) 
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regularization expressions like 



S=- 



l r 



( 7 " 7 " 7°^) 



(7) 



16 (7" 5) 



where P L R = (l±7 5 )/2. These 'counterterms' are just the unrenormalized operator 
(t$) inserted matrix elements, which form the basis for operator analysis that 
enters into RGE. In four dimensions the above expression in square brackets gives 
zero, but in n 7^ 4 dimensions, it has to be defined: is it a pole which must be 
subtracted, or is it a finite quantity which should not require a subtraction if one 
uses the MS scheme to renormalize the operators? This is the main difference 
between the naive scheme of Ref.[l] and the scheme of Refs.[2,3]. In Ref.[2,3], the 
authors looked upon such an expression as new structures and introduced evanescent 
operators for renormalizaton, while in Ref.[l] we did not. We shall come back to 
this shortly. 

Let us go back to Eq. (6) and rewrite it such that the operators are renormalized 
in the MS scheme. 



The l/m^ terms have been discarded. In this scheme too, we encounter expressions 
like those given in Eq. (7). Again, if one's intent is just to use the Zimmermann 
indentities as a technique to calcualte Fig. 4, then this is not a real problem: the 
expression — T r^ T, which contains evanescent structures, is found in 

the first and second terms on the right in Eq. (8) with opposite signs. This makes 
the sum of the two terms independent of the definition (i.e. prescription) of that 
evanescent expression. One has various choices to define the expression in Eq. (7), 
as long as one uses the same choice to calculate the first and second terms on the 
right in Eq. (8). Two possible prescriptions are: 




(8) 
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a) Formally, one can think of Eq. (7) as a pole that multiplies evanescent structures, 
and therefore treat them as counterterms for the operators (e.g. t$) which need 
renormalization. In this scheme, one introduces evanescent operators to carry out 
the renormalization of the operators, but some extra work is needed to obtain 
physical results. The renormalized Green functions of these evanescent operators 
must be included or be transformed away in the limit n — > 4. 

b) One can think of the expression in Eq. (7) as being finite and give a prescription 
to define its value. For instance, we can use the following definition: 

s = a (rp L ) 81S2 (rp L ) S3S4 (9) 

where a is arbitrary, but otherwise finite in the limit e — > 0. In this scheme, there 
is no need to introduce any evanescent operator. Only physical operators need be 
included in the analysis. We will present in section 4 a consistent scheme to evaluate 
the coefficient a. 

Let us summarize what we have said so far. The effective theory can be de- 
rived from the full theory with the use of the Zimmermann rearrangement: all the 
heavy effects in the full theory can be extracted by partitioning in momentum space 
the heavy graphs of the full theory. This is an alternative to the operator product 
expansion to construct the effective theory. In deriving the effective theory, one en- 
counters evanescent structures which must be defined in the context of dimensional 
regularization. This makes the renormalization of the physical operators artificially 
prescription dependent, and as a consequence the Wilson coefficients become pre- 
scription dependent as well. However, provided one stays being consistent in dealing 
with these evanescent structures, one is garanteed that £V CjT\^ ht (0[ en ) is inde- 
pendent of the prescription used to renormalize the various operators. One is then 
free to either include evanescent operators (assuming there are no inconsistencies 
within such a scheme) or not include them, in which case one must be ready to 
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present a consistent procedure in dealing with the evanescent structures. In sec- 
tion 5, we will treat in some details the processes b — > s in the LLA, and explain 
why one has to be careful when evanescent operators are first introduced and then 
transformed away by a finite renormalization. 

Before we close this section, let us complete this last example. For the process 
b + c — > c + s, in addition to Fig. 4 there are five more diagrams, plus external wave 
function renormalizations. If we follow the procedure to calculate the second term 
of Eq.(8) and its corresponding ones for the other diagrams, we obtain the effective 
Lagrangian to g 2 s 

LllL-c +s = -GcC 1 »P L s 8-fP L b 

+ <>*4?[™¥*™T* (io) 

- cR^P^b s^P^c Qlog^| - \ 

in which the evanescent operator is 

R K \1»P L ® ^R Xk P l = 7.7a7m^ ® 7 M 7 A 7 K ^ - 7*7a7^7 V ® YP L 

- ^P L ® 7„7A7^7 7" + 7m7a7^ <8> 7^7 7 M ^ , 

2 

G c = ^rV* s V c b, and the V's are the relevant Cabbibo-Kobayashi-Maskawa matrix 

w 

elements. This agrees with the result given by Buras and Weisz in Ref.[7]. The 
constant -11 in the first term of Eq.(10) is the one loop matching condition for the 
physical Wilson coefficient and it depends on the choice of the evanescent operator 
of Eq.(ll). It is needed for subleading logarithm approximation in solving RGE. 
Please note that the first term of Eq. (8) is also dependent on evanescent operators. 
It is only the S- matrix elements which should be and will be basis independent. 

We have also explicitly checked that the method is internally consistent up to 
two loop order in Feynman diagram evaluation of the processes b — > s"f and b — > sG. 
The results can be found in Ref. [4] . 
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3. Renormalization Group Equations (RGE) 



We want to give a quick review in this section, so that we can establish notation. 
We will define the renormalized operators and their anomalous dimension matrix. 
Then we will write down the RGE for the Wilson coefficients and that for the 
renormalized Green functions with the insertion of an operator. We will also discuss 
the finite transformations that take us from one renormalization scheme to another, 
and the advantage of using a mass independent renormalization scheme (MIRS) for 
the operators. The meaning of the leading logarithmic approximation (LLA) and 
the next-to-leading logarithmic approximation (NLLA) will be briefly reviewed in 
such a scheme. We will also say a few words about the matching of the effective 
theory with the full theory. It is understood that we use dimensional continuation 
as regularization. 

The renormalized operators are defined as linear combinations of a complete 
set of bare operators 



where the matrix Z is an infinite square matrix if we introduce evanescent operators, 
otherwise it is finite. It is determined order by order in perturbation theory by 
requiring that Y ren {O r l en ) be finite in the limit n — > 4. We would like to mention that 
the infinite matrices we will encounter have entries which eventually are unbounded 
(this will be discussed in section 5); we can therefore treat these matrices very 
formally. 

The /(/-independence of bare quantities allows us to derive a RGE for the Green 
functions with the insertion of a renormalized operator: 



bare 



(12) 




(13) 
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where denotes symbolically the anomalous dimension of the external fields, 

and 7 y is the anomalous dimension mixing matrix of the operators given by 

T« = -(^ Z «*) Z */- ( 14 ) 
The same mixing matrix also governs the evolution of the Wilson coefficients C l 

= CjTj*, (15) 

which is a result of the /(/-independence of CjOj. Eqs. (12-15) are valid for all 
schemes, including mass dependent renormalization schemes. They are also valid 
(at least formally) if we include evanescent operators in the analysis. 

A change in the renormalization of the operators must be compensated by a 
change in the Wilson coefficients. This amounts to making a finite transformation 
such that the sum of products CiO\ en is invariant. We shall drop the upperscript 
ren from now on. Under a finite renormalization (1+A), various quantities are 
transformed into the corresponding primed ones 

O' = (1 + A) O = (1 + A) ZO {bare) = z'O ibare \ 

C' = C(l + A)-\ (16) 

(1 + A)- 1 + (i + A) 7 (1 + A)" 1 , 



/ _ dZ> 
7 = 



d/u, 

where the matrix A is assumed finite but otherwise arbitrary, and may depend 
explicitly on fx. 

We will now discuss the advantages of using a mass independent scheme. In 
such a scheme, the Wilson coefficient vector C is assumed to have the following 
expansion: 

°° / 2 \ n 

C = zZW log ^ ) k 0) + ciV + • • • + + . . .1 (17) 

n=0 V ^ / L 
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where stands for the mass of the heavy field which is integrated out, and g = 
g(fjL). In Eq. (17) we have shown the explicit dependence of C on g, nih and fi; the 
c£f do not depend explicitly on g, ran and \i. In the LLA we keep only, in the 
NLLA we keep both c (0) and c (1) , etc ... 




To solve the RGE of Eq. (15) for C we need the initial (or boundary) conditions 
C{n)\ lJi = rnh . In a mass independent scheme, the initial conditions are given by 
C(m h ) = c ( q ] for the LLA, and C(m h ) = c|, 0) + g(m h ) 2 c { 1} for the NLLA. These are 
all very simple and clean in principle. As a contradistinction, in a mass dependent 
scheme where the have an explicit dependence on \i, the boundary conditions 
become C(m^) = J2^Lo 9( m h) 2p % { m h)i and one cannot truncate the series, say to 
the first term to obtain the LLA result. This is because each term g(mh) 2p Cq (rrih) 
contains "[sr 2 log(m);)] n " and therefore has to be taken into account to fulfill a 
consistent partial sum for the LLA. 

As we saw, in a mass-independent scheme for the NLLA, we need to com- 
pute Cq 1 '' as well to obtain the boundary conditions. We would like to repeat that 
these NLLA initial conditions (cq 1 ^) are scheme-dependent, i.e. they depend on the 
particular choice of a mass independent scheme one uses to renormalize the oper- 
ators. Therefore the C(ji) we obtain in the NLLA are scheme-dependent as well. 
Again, the sum Cj(//)Oj(/i) is scheme-independent. For instance, for four fermion 
processes, the NLLA results are given by: 

[cwoof " A = cr A WiY tree) + crnOiY 1 - 10 ^ . (19) 

This is illustrated in Fig. 5. In this regard, we differ from a remark in Ref.[8] 
by Dugan and Grinstein. In Eq. (19) the scheme-dependence of the (j NLLA w [\\ 
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be precisely canceled by that of r(Oi)^ 1oo p\ We have already given an explicit 
example at the end of the last section in this regard. 

In the renormalization group analysis of b — > s processes, all authors have 
assumed that their schemes are mass-independent. A closer examination seems 
warranted, however. 

4. Evanescent structures in the t'Hooft-Veltman scheme 

In this section, we present in the t'Hooft-Veltman scheme certain equations 
which are needed to define the renormalized graphs with an arbitrary number of 
external fermions, with or without insertions of operators. We are particularly 
interested in defining the expression in Eq. (7) without introducing any evanescent 
operator. 

In the t'Hooft-Veltman scheme, the external particles are kept in four dimen- 
sions: their momenta have the form p ext = (j>i,P2,P3, Pa, 0, . . . , 0), and the Lorentz 
indices attached to them run from 1 to 4. The spinors attached to these external 
particles will be elaborated. We will first solve the Dirac equation in even inte- 
ger n dimensions. Among the solutions we obtain, we will choose two to describe 
the two physical degrees of freedom of an external particle and two others for the 
anti-particle. These particular solutions will have the property that only the first 
four components of are non-zero. This will allow us to write for the spinors of the 
external particles as 

u(p)=Vu(p) , u(p)=u(p)V, (20) 

where V is some projector which will be defined in a moment. This projector, which 
was first introduced in Ref [9] , will be used to define expressions like those given in 
Eq. (7). 
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We start with n hermitian matrices 7 M (fi = 1, 2, . . . , n =even), satisfying 



(21) 



These are 2 n l 2 x 2 n / 2 matrices. Their construction is inductive and is explained 
in Ref. [10]. Let us briefly review the main ideas. The induction starts at m = 2 
(four dimensional Dirac matrices). We then assume that we have constructed 2m 
hermitian matrices 7^ 2m ); °f dimension 2 m x 2 m , which satisfy Eq. (21). This is 
extended to 2m + 2 hermitian matrices 7^ m+2 )' °f dimension 2 m+1 x 2 m+1 , which 
satisfy Eq. (21) as follows. We choose 



(A 4 
7(2m) 




7(2m) 



if 1 < jj, < 2m, 



and 



7 2m+l _ 
'(2m+2) 



where 



7(2m) 
7(2m) 



7(2m) 



2m+2 _ 
'(2m+2) 



il{2m) 
-*7(2m) 



•m l „,2m 
_Z 7(2m) • • -7(2m)- 



One can easily check that all the matrices 7^+2) are hermitian, and satisfy Eq. 
(21). 

From now on, we will write 7 M instead of 7^ . In this representation, the first 
four matrices 7*, i = 1, . . . , 4, have the special form: 

/ 7 W \ 



7 



\ 



(22) 



where 7^^ are the 4x4 Dirac 7-matrices, chosen such that 7^ is diagonal: 

Z 1 



4 

7(4) 



(23) 



V 
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The non-trivial solutions of the Dirac equation, i. e. those which do not vanish 
identically, which satisfy 

(i^.p + to) u(p) = 0, 

, , (24) 
p 2 + m 2 = 0, 

are easily found to be 



[VW] =6^ ; z,j = l,2,3,4,...,2-/ 2 . 



(25) 



Those with p 4 = +iE, E = \/p 2 + m 2 , will correspond to particles, while 
others with p 4 = —iE will be for anti-particles. For concreteness we will discuss 
the particle-solutions only, where the 2 n / 2_1 independent spinors have been afixed 
with a normalizing factor, 

"" >(p) = 71ICT ( - iTP + m)yW ' (26) 
1=1,2, 5, 6, 9, 10, . . . , 4£ + 1, U + 2, . . . 

< £ < 2 n / 2 " 2 - 1. The normalization of the spinors is: 

(w (<) ) f = 1, for eachi (27) 

and one can easily derive: 

U W = for each i, 



2£ 



where z takes the values 1, 2, 5, 6, . . . , Ai + 1, 4£ + 2, . . .. In deriving the above equa- 
tion, we have used 

2 

i 

where the upperscript * stands for transpose. 
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If we go to the rest frame of the particle (p = and p 4 = +im) , the spinors in 
Eq. (26) become: 

u {1) = (1,0,0,0,0,0,0,0,0,...,0), 
u {2) = (0,1,0,0,0,0,0,0,0,...,0), 

= (0, 0, 0, 0, 1, 0, 0, 0, 0, . . . , 0), (29) 

= (0,0,0,0,0,1,0,0,0,...,0), 

In the real world where n = 4, a spin-i particle has only two degrees of freedom. 
This nomenclature is kept in n dimensions where we choose to designate the spinors 

and given in Eq. (29) as the two physical degrees of freedom of a particle 
at rest: for a particle with spin up, and with spin down. The other spinors 
will describe non-physical degrees of freedom which must disappear when the limit 
n — > 4 is taken. Similarly, we will choose u^ 3 \p) and u^ 4 \p), p = 0, p^ = —im, to 
describe the anti-particle at rest, with spin up and down, respectively. 

Therefore, when computing Feynman diagrams in n dimensions, because the 
external fermions are kept in four dimensions, they will be described by states (or 
spinors) which are linear combinations of u^(p) and u^ 2 \p) . Such states, which 
we will call physical spinors, have the form: 

u(p)p hys = a u {1) (p) + 13 u {i) (p) = (x, x, x, x, 0, ... , 0) (30) 

i.e. only the first four components of u p h ys can be non-zero. This can easily be 
checked if we note that when p is the momentum of a particle in four dimensions, 
then the matrix 'j.p is of the form 

4 / 7(4) -P 



1-V = J2 7 M P" = 



7(4) -V (31) 



V 



and therefore the vectors {—i^.p + m)V 1,2 are of the form shown in Eq. (30). We 
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now define a 2 n / 2 x 2 n / 2 matrix V 



( n n)- (32) 





In the above equation 7 4 is the 4x4 identity matrix. V is a projector (VV = V) 
which satisfies: 

U(p)phys = V U(p)phys , u(p)phys = u(p)phys V . (33) 

We can now state our prescription in handling gamma matrices: "when com- 
puting Feynman diagrams in the t'Hooft-Veltman scheme, we simply attach a factor 
V to each external fermion" (see Eq. (32)). For instance, for four quark processes, 
Eq. (7) will be rewritten as: 

S = -\ ( V r l v l a P L V) si S2 ( V ^ Y 7 a P, V) S3 S4 

(34) 

-i6(vrp L v) slS2 (vrp L v) S3S4 \- 

We will see later that the presence of V will make the above expression collapse 
into Eq. (9), in fact with a=0. In other words, the gamma algebra in the first term 
may as well have been done in four dimensions. 

One can check that for n even, V can be written as: 

n/2 x (35) 

p=3 ^ ^ 

This is easily seen if we observe that the products of matrices 72 P -i72 P , P > 3 are 
diagonal; if I m and m denote identity and the null 2 m x 2 m matrices, respectively, 
we have: 

1 ~ hn-lln _ ( / 2 ™/2-l 



(36) 
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Using Eq. (35), we easily derive: 

[7>,7 M ]=0 if n= 1,2,3,4, 

VYV = if ^>4, (37) 
[P,7*] =0. 

The last identity in Eq. (37) holds if 7- anticommutes with all 7 M . It also holds 
in the t'Hooft- Veltman definition, where 7- anticommutes with the first four 7- 
matrices but commutes with the remaining n — 4 matrices. This matrix should not 
be confused with the matrix 7 M , when /j, = 5. Moreover, whether 7- is defined as 
the product of the first four 7-matrices, or as the product of the n 7-matrices (n 
even), we have: 

f P-l P-2 P-3 P-4 

V^V = V 7V7V V = - V 7 ^i 7 ^ 7 ^3 7 M4 ( 38 ) 

where the barred indices take on the values 1,2,3 and 4 only. 

As far as four fermion processes are concerned, we need to evaluate the following 
quantities: 

Tij^iVSiP).^ (VSiV) asat (39) 

where the Si are products of 7-matrices, and possibly 7- s. The algorithm to com- 
pute such expressions is the following: we first use the sum splitting notation 

(7% S2 (7a 3S4 = M S1S2 (7 A ) S1S2 + (7 A ) S1S2 (7 A ) S1S2 , (40) 

where 

H = l,...,n, 

^ = 1,...,4, (41) 
(x = 5, . . . , n. 

Then we move all 7^, to the left with the use of Eq. (21). If 7^'s with the t'Hooft- 
Veltman prescription are present, they are also moved to the left. [For that matter, 
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one could use a completely anti-commuting 7- if there are no traces involving 7-] . 
At this point, %j has the form 

%j 7 W 2 • • • Y m (?VV 2 . . . T*>V)] ai S2 

(42) 

[ ^ . . . 7 5 7 a l7 a a 7 a, . . . ^p)] ^ ^ . 

In the expression above, we do not have any hatted index which is not summed, 
because all the external particles, and hence free indices, are in four dimensions. 
The following identity is then used to reduce Eq. (42): 



S1S2 \ / S3 s 4 

(-l)W2 m! ( w/2 -2)! (43) 

- W» - (n/2-2-m/2)!(m/2)! lf ™ ^ 
if m odd. 

where the square brackets denote antisymmetrization in all m enclosed indices (in- 
cluding the usual 1/m! factor). At this stage, the remaining Dirac 7-algebra is 
performed in four dimensions to obtain %j. 

Let us be reminded that Eq. (43) has been obtained for even integer, and is 
used as a definition for complex values of n. It is easily obtained if we observe that: 

V H72 P -i72 P ) V = V for p > 3 (44) 

As an example, an application of Eq. (43) gives 

{ V ^ ll ^ 2V ) SiS2 (?W 2 ^) S3S4 = (45) 

where there is no antisymmetrization over i>\ and z> 2 . The final result for %j must 
then have the form: 

5 

%i ~ E *ab (rj slS2 (r B ) S3S4 , (46) 

A,B=1 

where 

T A = V, V^V, Va^V, V^l-V, V^V (47) 
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This is obvious because any 7-string VSV can be written as 

vsp= ( S o )' (48) 

where S4 is the projected 4x4 submatrix. Thereupon, it can be written as a linear 
combination of the r a listed above. One can also perform " four dimensional Fierz 
transformation" for the T A 's if that is necessary. 

Before we close this section, we would like to make a few more remarks. 

1. The final result for T is not fully Lorentz invariant (although it is invariant under 
a " four dimensional" Lorentz transformation) . One could however rewrite Eq. (46) 
with the use the following relations: 

(49) 

Va» v V <g> Va^V = Va^V <g> Va^V + a V <g> V 
where ciq is proportional to e. This apparently "restores" the full Lorentz invariance 

(apart from the projector V which can then be reabsorbed into the external spinors: 

Vu(p) = u(p)). 

2. The spinors u^'(p) and uS 2 \p) which describe external particles, satisfy: 

2 

]T M (0 (p) u W (p) = 1 (-z 7 .p + m) pl±^ip (-i 7 .p + m) 

2£(m + £) 2 (5Q) 

= P ^ (-i 7 .p + m) T 3 
The above relation is easily obtained if we observe that p = (pi,P2,P3,P4:, 0, . . . , 0), 
and 

J^V^V^ 1 = V^-^V (51) 
Eq. (50) is used to calculate cross-sections, as an example. 

3. If we have more than four external fermions, we have to generalize Eq. (43). For 
example, expressions like 

(v 7 [Al . . . 7 A2m] V} <S> (v 7^ . . . 7 i>2 ™ 1 ® (v 7 [Al . . . Y 2m] V^j (52) 
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are first calculated in even integer n dimensions. The result is then taken as the 
definition of that expression for complex values of n. The Lorentz indices in the 
above equation must of course come in pair. Here are a few examples: 

^-p^IAi^^As^A-i] pj (g, ^p^^^As^] pj ^p 7 [£1^2 ^1^2] pj _ 

-3(n-4)(n-6)(n-8) V®V®V, 

^p 7 [Ai 7 A2 7 A3 7 A 4 ] pj jp 7 [Ai 7 A 2 ] pj ^p 7 [A3 7 ^4] pj = (53) 
(n-4)(n-6) V®V®V, 

^p 7 [Ai 7 A 2 ] pj (g, ^p 7 [A2 7 A 3 ] pj (g, jp 7 [A3 7 Ai] pj _ 0. 

4. We now give an example at the two-loop order to show how the projector V 
can be used to derive the effective theory to this order. This is illustrated in Fig. 
6. Let T represent the two-loop diagram on the left , 70 its tree (heavy) subgraph, 
and 71 its one-loop (heavy) subgraph. 70 is represented by the diagram on the left 
in Fig. 3, and 71 by the diagram on the left in Fig. 4. Since T is finite, the Dirac 
7-algebra can be done in four dimensions with no ambiguities. We can also attach 
a factor V to each external fermion without changing the value of the graph in the 
limit n — > 4. The Zimmermann algebraic rearrangement for F gives: 



r = 



1 



1 - T W) r(°) 

71 / 70 



+ 
+ 
+ 



M 0) ) 



-(0) 

71 



7", 



(0) 



(54) 



M 0) ) (l-T«>)(l-r(?) 



r 



A simple power counting argument shows that the last term in the above equation 
is of order and will be discarded. The first, second and third terms in Eq. 

(54) can be identified to the first, second and third diagrams on the right in Fig. 6, 
respectively. Note that Eq. (54) uses subtraction at zero momentum to renormalize 
the four-quark operators. Also, each of the first three terms are finite: the Dirac 
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algebra can be done in four dimensions, and the external fermions carry implicitly 
a factor V . Let us now rewrite Eq. (54) such that the four-quark operators are 
renormalized in the MS scheme: 



What we have done is simply to rearrange the terms in Eq. (54). In the above 
equation, " lim e ^o [• • •] " means precisely that we compute the finite subgraph " 
[. . .] " by first assuming its external particles are in four dimensions. The result is 
then reinserted in n dimensions. The operation extracts the pole of the graph 
(or subgraph) 7 in the same fashion. Let us recall that external particles carry 
a factor V: when the result is reinserted in n dimensions, the factor V must be 
removed. 

When this projector method is used to renormalize the various four-quark op- 
erators, together with the use of an anti-commuting 7 5 , the LLA results of b — > s 
processes are identical to those of Ref. [1] . 

5. b — > s processes in the LLA 

In this section we discuss in the context of dimensional regularization and 
evanescent operators the processes b — > s, in the LLA. It is independent of section 
4, and no projector is used. Here we present in the MS scheme a two loop calcula- 
tion of the mixing of all the evanescent operators with the physical dimension five 
operators (also called the dipole moment operators). The results in other schemes 
can be obtained with the use of Eq. (16). We want to take note that we are now 
dealing with an infinite set of operators. It follows that the dimension of the mix- 
ing matrix is infinite. In order to make mathematical sense, in principle one needs 





(55) 
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to address the issues of boundedness, convergence, asymptotic completeness, etc. 
These unfortunately we do not know how to deal with because some of the matrix 
elements are in fact unbounded in value. An example of inconsistency is given to 
show the potential pitfalls. On a formal basis, however, we shall show explicitly that 
there is a similarity transformation which will decouple the evanescent set from the 
physical operators we are interested in. Thus, based on this 'physical argument', 
one may give meaning to the mixing matrix as such. 

For b — > s processes, the commonly used basis of physical operators is 



Oei = 






62 = 


(s a rP L cp)(c^P L b a ), 




063 = 


(s a YP L u a ) (u^P L b p ) 7 




o m = 


(s a ^P L up) (up^P L b a ), 




o G5 = 


{s al »P L b a ){u f3l »P L U l3 + . 


, + b^PM^ 


o m = 


(s a ^P L bp) (u/3^P L U a + . 


. + b p ^PM, 


o G7 = 


(s a ^P L b a )(up^P R Uf3 + . 


. + 1^pm, 


o 68 = 


(s a ^P L bp) {U{3^P R U a + . 


. + b p ^p R b a ), 


051 = 


^ig s sG^a^ u (m s P L + m b P R )b, 


52 = 


]^ieQ d sF tlv a^{m s P L +m h P R )b ; Q d 



a and (3 are color indices. 

A few words about this choice of basis are in order. When the operators 
(965, 068 are inserted in the Green functions corresponding to the processes 
b — > S7, sG, one encounters traces of Dirac matrices and 7 5 . This has led some 
authors to use a scheme different from the naive dimensional regularization scheme 
(i.e. 7 5 anti-commuting), such as the t'Hooft-Veltman scheme or the dimensional 
reduction scheme. It should be observed that when we include evanescent operators, 
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this 7 5 problem is artificial and is created by the choice of the basis listed in Eq. 
(56). One can in fact choose a different basis of physical operators and use an 
anticommuting 7 5 unambiguously. For instance, one can replace the operators 
Oq^ , . . . , by the following new operators 

<?65 = {s a ^P L b a ) (UfB^Up + ... + IfB^bp), 
<?66 = {San^PM {u/3^U a + ... + l^b a ), 

(57). 

Q 6 7 = ( WW'^a) ( WtVW + • • • + 6/37 7^7 

Q 6 8 = {s al [ ^r ] PM (w^VW + ■■■ + WW&a), 

which do not generate any trace involving 7 5 . Thus if we include the corresponding 
evanescent four-quark operators (see below), then the use of a 7 5 anticommuting 
with all the 7 M is justified. 

Following the notation of Ref. [8] , we define the totally antisymmetric matrices 

ryW] = 7 lMl ... 7 A»p] j (58) 

with the understanding that is the identity matrix. The basis of evanescent 
operators that can mix with (567,68 can be chosen as (below p is odd, and p > 5): 

Qll = (Sal [ ^ ] PM (%7 W «/3 + • • • + &/37 M M, 

(59) 

^68 = (Sal [ ^ ] PM (*W MP V + • • • + h-f^K). 

We have checked explicitly that these are the only evanescent operators that can 
mix with the dimension five operators O5152, in the MS scheme. 
Let us first rewrite Eq. (12) as: 

ren = J- Z t3l i- eD > 0) are (60) 
j 

where the matrix Z is dimensionless, D l = if 1 corresponds to 051,52, and D l = 1 
if 1 corresponds to a four quark operator (physical or evanescent). The factor fi~ eDl 
is introduced to make the mass dimension of 0[ en constant: dim (O z ) = n + 2 for 
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all i. With this choice, all the Wilson coefficients have dimension mass 2 . To order 
g 2 , we have in the MS scheme: 



1 + 



9 2 s 



8tt 2 



■B, 



7* J = 



8n 2 



(61) 



(1 + D t -D 3 )B. 



We want to calculate the matrix B. As discussed, because we will not encounter 
any traces involving 7 5 , we will take for definiteness 7 5 anti-commuting. Also, the 
calculation will be done in the general linear covariant gauges " ^(d^G^) 2 " for 
the gluons. 



First, we insert the evanescent operators in a four quark process: this gives the 
mixing of these operators with themselves. There are 6 one loop diagrams which 
give: 

1 32 8 

B Qlr ■ Qlr =6 6p+2 > q + (8 " Y P + 3 6p > q 

+ (~5p+ y P 2 - 2p 3 + 5 P _ 2)(J , 
1 41 1 

Bqi & , Q%7 = ~ \ S P+ 2, g + (3 - 6p + ^p 2 ) S p , q (62) 
+ ( y P - y P 2 + 3p 3 - - p 4 ) 6 p -2, q , 

.35 287 2 ^ 3 7 4 . „ 
+ ( Y P-^P 2 + 7p*--p 4 )5 p _ 2 , q . 



The above results will be needed when we insert the evanescent operators in b — > S7 
and 6 — > sG, to order g 2 . They have been obtained with the help of the following 
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identities: 

7 « 7 K1 7 « = (-l)P(n - 2p)^\ 

ry°<^lVp] ^ 7 « 7 K] _ ^p]^ ^ -yM^ _ 

7 K+i] ^ 7 K+i] + p ( n + 1 _ p ) 7 K-i] ^ 7 [Mp-i] ; (gg) 
7 a 7 [^p] 7 [A*p] 7 a = 

(-l) p ^>p+d ® 7 K+i] _ p( n + l - p) 7 K-il 7 K-i] j 5 

which were also derived in Ref. [8] . 

Next, we insert Qq 7 68 in b — > S7 and 6 — > sG to the lowest order. This is 
illustrated in Fig. 7. The results are finite, and will be used (at the end of the 
calculation) to " renormalize away" these evanescent operators; in other words, we 
will redefine them such that they vanish at the lowest order (for the LLA). For 
this, we define a finite matrix A such that if Oi is evanescent, then 

r(O t + A tj O J ) b ^ sl ' sG ^0 when e -> 0. (64) 

The above equation defines A at lowest order in the strong coupling constant, and 
Oj = 051,52- 

Before we present the results, we would like to give the various identities which 
are used to simplify the n-dimensional Dirac 7-algebra. We can easily derive, for 
arbitrary p, 

MM = (_!)p(p-i)/2 n! (65 ) 

[n — p)\ 

For p > 5 and n — > 4, Eq. (65) becomes 

7 M 7 [mp] = _ 24e (_i)p(p+i)/2 (p _ 5)! + Q ( £ 2) _ ( 66 ) 

More generally, we have (see also Ref. [11]): 

7 K1 7 K1 7 K] = 7 K ¥ Mp] 0(p) k] x 67 ) 

' ' ' 7 [«q] 7 [«q] 3 ' V ' 
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where the parameters j3q are given by the following recursion relation: 



4 P) =1 , P[ P) = (-l) p (n-2p), 
«Cp) = (_i)H^-i (n _ 2p) M + {q _ 1)(n _ q + 2) ^(p) 



(68) 

V/-2 



For the present calculation, the values of for q = 0, 1, . . . , 5 only, are needed. 
The above recursion relation is obtained with the help of the following identity: 



(69) 



4 0) = = 1 , a?> = 0, 
a[ p) =0 if r < or if r > p, 
aOH- 1 ) = p( n + 1 - p)^- 1 ) + (-lfa^V 
Also, the following relation holds: 

/#°=£ 4 a) [(-l) p (n-2p)] r . (70) 

r=0 



Our results for the matrix A are, when p > 5: 



A n p r>_ = ^4r>p r>„ = — A f 



Hp) 



(71) 



A Q P 68 ,0 51 = °, 

where 

h(p) = (_1)pCp+i)/2 (p _ 5)| (p _ !)(p _ 3 ). ( 7 2) 

Next, we insert these evanescent operators in 6 — > S7 and 6 — > sG to order 
<7g. This is illustrated in Figs. 8, 9 and 10. To extract the infinities of the two 
loop diagrams, we need to calculate the double pole (1/e 2 ) only; this is because the 
Dirac 7-algebra always generates an extra e, as can be seen from Eq. (66). In the 
MS scheme, it is a well known fact that the double pole of a two-loop diagram is 
a polynomial in the external momenta and masses. One can therefore expand each 
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diagram in a power series in the external momenta and masses, to order mass 2 . The 
calculation is straightforward and is done with the help of the symbolic manipulation 
program Schoonschip Ref.[12]. The results for 7^ (see Eq. (72)) are: 

9 2 s w w 40 28 7 2 
7q£ 7 ,o 51 = ^4 Hp) (—3- + yP" ), 



9 2 



(73) 



^o 5 , = ^Hp) (-2 + 8p-2^) 

^ 40 64 16 2 , 

7Q- 7 ,0 52 = ^4 h(p) (-y + yP-yp) 

7Q- 8 ,o 52 = ^(-8)h(p). 

Let us imagine that we solve Eq. (13), which formally has the solution 

CO*) = C{m w ) T exp ( f ^ 7 (/i) ) , (74) 

where T stands for normal ordering of the matrices 7 (^). For the LLA results, the 
mixing matrices 7 (^) commute for different fx and therefore the T-product can be 
omitted. Thus We obtain: 

C(fi) = C(m w ) exp ^ ^L- log(^) ^ , (75) 



where 

7= 8^ 7 ' 



>lw — 9/ \ — 1 1" 1 „ 9 y s 



(76) 

(^) log 



g 2 (m w ) " ' 16tt 2 ^ svrv ~° ^ 2 
In the above equation, the only non-vanishing boundary conditions needed in the 

LLA are C 0eii63 (m w ) and Co 51i52 (m w ) (see Eq. (56)). 

The matrix 7 has entries that become large as p becomes large (see Eqs. (62 
and 73) . This tells us that the evanescent Wilson coefficients may become large as 
p — > 00, and the convergence of 

^2 C Oi10i,0 51 , 



i=evan 
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is immediately brought into question. For example, we do not know how to evaluate 
numerically Cq£ 7 at say fx = 5 GeV. 

The handling of such an infinite set of operators is quite delicate. To illustrate, 
suppose one chooses instead a new basis {Ui} of evanescent operators, defined by: 

Uqt = Q\ 7 , U 68 = Qg 8 , 

U & 2 = + (p+1 y p . 4) Qif fol ' P ^ 5 - (77) 
- «. + ( p+ l) 1 ( p -4) «&' for „ > 5. 
This new basis has the property that in the limit e — > 0, the renormalized Green 
functions of U% 7 68 , p > 7, vanish at the lowest order for 6 — > S7 and b — > sG. 
Only Uq 7 68 have a non-vanishing value. Because the basis {t/i} is infinite, one may 
querry whether the operators Uq 7 68 are in fact linearly independent of the other 
Uq 7 68 , p>7. One may be flippant and formally writes 

00 

U$> = J2d P U% 7 +2 (78) 
where the coefficients d p are obtained by the following recursion relation: 



1 

(p+lXp-4) 



<^p+2 = -7 — — tt7 tt rf P ; d 5 = 1. (79) 



Is one allowed then to choose a reduced basis {Ui} which does not contain the two 
evanescent operators Uq 7 68 ? Formally, it seems that nothing should prevent one 
to start with this reduced new basis to renormalize all the physical and evanescent 
operators. The question is then: do the full basis and the reduced basis give the 
same on-shell results? If one uses Eq. (78), then the Green functions to b — > S7 of 
Uq 7 at the lowest order are not well defined, since the term on the left in Eq. (78) 
gives non-zero, and each term on the right gives formally zero. This "problem" is 
obviously related to the fact that one has an infinite number of evanescent operators. 
It would not exist if the number of evanescent operators were finite. Unfortunately, 
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as we have seen, one cannot carry out the program of renormalization with a finite 
number of evanescent operators, except for a prescription in which no evanescent 
operator is introduced altogether. 

One may skirt this potential difficulty by the following: the LLA result for the 
amplitude of b — > S7 is 

[C i r(O l )]" A = C- LA T(O t ) lowest order (80) 

where the sum is over all physical and evanescent operators. If we use Eq. (64), 
then we can rewrite Eq. (80) as 



LLA ^ — , 

[c,r(0i)] = E 

i=phys 



.^lowest order (81) 



j=evan 

The above result was obtained in the MS scheme. It is however the same in all 
schemes (including mass-dependent schemes). One has redefined the evanescent 
operators such that their physical matrix elements vanish at the lowest order. This 
leads to 

Z -> Z' = l + A+-^-B', 

8 7r 2 e ' (82) 

= (1 + A) Z. 

This manipulation was implicitly used in Ref. [2,3] and the anomalous dimension 
matrix becomes (see Eq. (16)): 

7-7'=(l + ^)7(l-^), (83) 

where 7 is the mixing matrix obtained in the MS scheme. To the order, we have 
used (1 + A) -1 = 1 — A , A 2 = 0. One can easily check that 



(84) 



68^52 
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Effectively, the mixing matrix between the four quark and the dipole moment op- 
erators is collapsed into a finite dimensional one, which has elements from physical 
to physical. 

We can now finish constructing the mixing matrix between the physical four 
quark operators and the dipole moment operators. They agree with the results 
by Ciuchini et. al. Parenthetically, we remark again that in view of our previous 
discussion this requires the assumption that lim e ^o J2i= e vanes CiT(Oi) = 0. In our 
opinion, this has not been firmly established because, as we have seen, Cqp^ may 
become unbounded as p approaches infinity. 

The renormalized physical operators in the above scheme (where evanescent 
operators are first introduced then renormalized away) are related to those in the 
naive scheme in which there are no such operators by a finite transformation which 
depends explicitly on fx: 

[O phys r- scheme = (1 + g 2 s log(p/m H9ht ) F + . . .) [0 phy8 ] naive sc/ ™, (85) 

where the finite matrix Fij ~ "fi^Akj. This means that the Wilson coefficients 
cannot be of the form given by Eq.(17) in both schemes. One of the schemes must 
be mass-dependent. This is what we referred to at the beginning of this article as the 
light-mass dependence of the Cj. Thus, the problem has been shifted to a problem 
of factorization, namely, finding a prescription in which the Wilson coefficients do 
not depend on the masses of the light particles. This is our present understanding 
to account for the difference in results between Ref.[l] and Ref.[3]. One has to go 
to at least three loop order to settle the issue. 

6. Concluding Remarks 

In this article, we have attempted to locate the cause for the difference in 
the short distance analysis of b — > s'j between the scheme in Ref.[l,4] in which no 
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evanescent operators are introduced and those schemes in Ref.[2,3], which contain 
such operators. We have argued that if one follows a systematic formulation to 
arrive at an effective Lagrangian, which for our preference is the Zimmermann 
rearrangement, then there is ample freedom in divising different schemes. We have 
explained how the Wilson coefficients Ci depend on the scheme used to renormalize 
the operators. The boundary conditions also depend on the procedure. This can 
be easily understood if we observe that the identification and the subtraction of the 
"infinities" are done differently in different schemes. Consequently, the Ci to that 
order must be different, which affect the matching conditions. Such difference will 
exactly cancel out the difference in the anomalous mixing between various schemes. 
It is the S-matrix elements of the product Yl CiOi which is independent of the 
scheme used. We have also argued that the problem of 7 s is artificial, since all 
quantities are finite to begin with; provided one is consistent all prescriptions for 
handling 7 5 are equivalent. 

In section 4, we have elaborated on a version of the t'Hooft Veltman scheme, in 
which all external particles are kept in four dimensions and we have developed an 
algorithm to handle the Dirac algebra for processes with more than two fermions, in 
particular when four-quark operators are introduced into the analysis. This scheme 
as laid out so far may, however, be incomplete, since it breaks Lorentz invariance in 
n dimensions. It would be instructive to construct a Lagrangian which is Lorentz 
invariant only in four dimensions; for instance, we could replace every fermion field 
ip by Vip (and ip by ipV) at the lagrangian level . This would have the advantage 
that in the evaluation of Feynman graphs derived from such a Lagrangian, the Dirac 
7-algebra could be handled in four dimensions. We do not know if such a theory 
(or scheme) is consistent, but from a theoretical point of view we believe it is worth 
investigating. 

In section 5, we have followed the approach in Ref.[2,3] in introducing evanes- 
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cent operators. We have shown that some entries in the mixing matrix are un- 
bounded. However, by a formal finite transformation, which is induced by impos- 
ing a requirement that the matrix elements of the evanescent operators vanish at 
the lowest order, the mixing matrix elements of the new evanescent basis with the 
magnetic operators formally vanish also. The results of Ref.[3] (as far as the anoma- 
lous dimension mixing matrix of the physical operators) have been reproduced. We 
nonetheless are concerned that this scheme lacks certain justification in some for- 
mal manipulations which are related to the convergence of the infinite series with 
coefficients that eventually go to infinity. 

We have also stressed that in going from finite order analysis to LLA with 
the help of RGE, we need to make sure that mass independence of Ci is valid to 
facilitate the summation of leading logarithms. By this, we mean that the only 
log// dependence in Cj can only be in the form log(M/ ieauy / / u). We have not been 
able to establish this property in either scheme. We however believe that this is 
the most important aspect which needs to be clarified: it is our contention that 
this should account for the difference between the results of Ref. [1] and Ref. [3] and 
an unambiguous way to settle this is to perform an explicit higher loop Feynman 
diagram calculation, where the discrepancies first show up. Further work needs to 
be done. 
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FIGURE CAPTION 

Figure 1. Diagrams contributing to b — > S7 at lowest order in the Standard Model. 

Figure 2. The Zimmermann rearrangement for a one- loop graph contributing to 
b — > S7. 

Figure 3. The Zimmermann rearrangement for four-quark processes at the tree level. 
The effective vertex on the right is reinserted in one-loop diagrams contributing to 
b — > S7, as shown on the first graph on the right in Fig. 2. 

Figure 4. An example of the Zimmermann rearrangement for four-quark processes 
to order g\g\. 

Figure 5. Graphical representation of the next-to- leading logarithmic approximation 
final results for four-quark processes. In this approximation we are summing all 
terms of the forms g 2 s n log n (m^,//x 2 ) and <7 2n+2 log n (m^/// 2 ). The final results are 
prescription independent since the "ambiguities" (concerning 7 5 and the inclusion 
or exclusion of evanescent operators) in the graphs on the right exactly cancel. 

Figure 6. Example of the Zimmermann rearrangement to order g^g^ for be — > sc. 
In the LLA, one keeps the third graph on the right only; in the NLLA, one keeps 
the second graph on the right as well; in the next-to-next-to-leading logarithmic 
approximation, one must also include the contribution of the first graph on the 
right. Again, all the ambiguities we encounter in the graphs on the right cancel 
exactly, since the graph on the left is finite and prescription independent. 

Figure 7. Contribution of four-quark operators to b — > S7, sG at lowest order. 
These diagrams also give the finite mixing matrix A which is used to "renormalize 
away" the evanescent operators. 

Figure 8. Two-loop graphs contributing to the mixing of four-quark operators with 
the dipole moment operators. The external wavy line represents a photon or a 
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gluon, and the black square represents an insertion of a four-quark operator. 

Figure 9. More two-loop graphs contributing to b — > sG with an insertion of a 
four-quark operator. 

Figure 10. Graphs needed to remove the sub divergences of the two-loop graphs in 
Figs. 8,9. The crosses represent the necessary one-loop counterterms. 
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